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Abstract. Motivated by the PAMELA anomaly in the fluxes of cosmic -ray e+ and e^, we study 
the cosmic 7-ray induced by the inverse Compton (IC) scattering process in unstable dark matter 
' (— I , scenario assuming that the anomaly is due to the emission by the decay of dark matter. We cal- 

Q^' culate the fluxes of IC-induced 7-ray produced in our Galaxy and that from cosmological distance, 

and show that both of them are significant. We discuss a possibility that large dark matter mass over 
TeV scale might be constrained by the 7-ray observation by Fermi Gamma-ray Space Telescope. 
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^ : INTRODUCTION 
00 

^ ■ Recent observations of the fluxes of high-energy cosmic rays have made an impact on the 

^ ! understanding of the nature of dark matter (DM). In particular, the PAMELA experiment 

has observed an increasing behavior of the positron fraction in the cosmic ray in the 
^ ■ energy range of 10 GeV < Eg < 100 GeV (with Eg being the energy of e^) [1], which 

cannot be explained if we consider the conventional fluxes in astrophysics. This fact 
suggests that there may exist a non-standard source of energetic positron (and electron) 
^ ■ in our Galaxy. 

^ . One of the possibilities is unstable dark matter. If dark matter has lifetime of 0(10^^ — 

10^^ sec), and also if positron is produced by the decay jthe PAMELA anornaly may be 
explained; for early attempts, see, for example, [0,lli,ll,Sl7L[8Li,[il[nl[il[il 
14ll . (Other possibilities of explainin g th e PAMELA anomaly include the dark- matter 
annihilation [H HI 113, S [S I^IIM B and the positron emission from 



pulsars (2^.) In addition, a precise measurement of the total (e+ + flux has been 
performed by the Fermi Gamma-ray Space Telescope [26], whose results suggest that 
the flux of (e+ + ^ ) is proportional to ~ E^^. Although the interpretation of the Fermi 
result is still controversial, it has been discussed that the observed spectrum may be 
too hard to be consistent with the prediction of conventional astrophysical model, and 
that the e"^ observed by the Fermi experiment may be significantly contaminated by 
the produced by the decay of dark matter. Such a scenario suggests the mass of 
moM ~ 0(1 TeV) and the lifetime of Tdm ~ 0{lO^^ sec) S S IM • 



If the decay of dark matter is the source of the extra positrons observed by the 
PAMELA experiment, the emitted positron and electron produce photon via synchrotron 
radiation and inverse Compton (IC) scattering. In our Galaxy, energy loss rates of the 
energetic via these processes are of the same order, but the typical energy of the 
photon emitted by these processes is different. Since the magnetic field in our Galaxy is 
expected to be 0{l /iG), the energy of the synchrotron radiation is typically 10^^ eV 
when the energy of is 0(1 TeV). (The synchrotron radiation from the Galactic 
center is discussed in i30l fsill.') On the contrary, the IC process produces y-rays with 
higher energy. If an energetic with E^r^ 0{ \ TeV) scatters off the cosmic microwave 
background (CMB) photon, y-ray with ~ 0(1 — 10 GeV) is produced. In addition, in 
our Galaxy, there exist background photons from stars, which have higher energy than 
the CMB radiation. The IC scattering with those photons produces y-ray with higher 
energy. Importantly, the high-energy y-ray flux can be precisely measured by the Fermi 
telescope. Thus, in order to examine the scenario in which the PAMELA anomaly is 
explained by the decay of dark matter, it is important to study the energetic y-ray emitted 
by the IC process. 

In this paper, we study the flux of y-ray produced by the IC process in the decay- 
ing dark matter scenario. We pay particular attention to the parameter space in which 
the positron fraction is in agreement with the PAMELA results. In our previous work 



[1320, IC-induced y-ray in various direction from Galactic and extra-Galactic region are 
shown; we have seen that the extra- Galactic contribution can be as large as Galactic 
one. In particular, if we see y-ray in the direction off the Galactic center, signal of the 
IC-induced y-ray from extra-Galactic region may be observed. Here, paying attention 
to the observation data by Fermi, we compare theoretically-calculated y-ray with the 
observations. It will be shown that the parameter region moM ^ 4 TeV, which is sug- 
gested by PAMELA and Fermi observations, might be constrained by the observation of 



isotropic y-ray by Fermi [|33l1 . We will also show that the flux of the IC-induced y-ray 
in our Galaxy may be less than the observed flux in intermediate Galactic latitude (IGL) 
given by Fermi [|34ll when Tdm ~ 0(10^^ sec) to explain the PAMELA anomaly. 

FORMULA OF 7-RAY IN INVERSE- COMPTON PROCESS 

We first discuss the procedure to calculate the y-ray flux. The following formula is based 
onS. 



The total y-ray flux is given by the sum of two contributions: 

4.^^=4>(p^'-y)+4>S,^-'^°\ (1) 

where the first and second terms in the right-hand side are fluxes of y-ray produced 
in our Galaxy and that from cosmological distance (i.e., extra- Galactic contribution), 

respectively. Notice that is isotropic, while i^'^^^^^^^^ depends on direction we 

observe. 

In order to discuss the IC-induced y-ray, it is necessary to understand the spectrum of 
the parent e^. In our Galaxy, energetic is approximately in a random-walk motion 
because of the entangled magnetic field. Then, the energy spectrum fe (i.e., number 



density of (e+ + e ) per unit energy) in our Galaxy is described by the following 
diffusion equation: 

Ke{Ee)V^fe{E,x) + ^ [bioss{Ee,x)fe{Ee,x)] + Qe{Ee,x) = 0, (2) 

where Ke{Ee) is the diffusion coefficient, b\o^s,{Ee,x) is the energy loss rate, and 
Qe{Ee,x) is the source term. In considering the long-lived dark matter, the source 
term is given by 

"^DM ''^DM dt^e 

where p^^^^^^ is energy density of dark matter and dNg/dEg is energy distribution of 
from the decay of single dark matter. In our study, we adopt the isothermal halo 



density profile Il35h 



P^^^'^'\r)=PrJj^. (4) 
'core I ' 

where pQ ~ 0.43 GeV/cm-^ is the local halo density, rcore — 2.8 kpc is the core radius, 
tq ~ 8.5 kpc is the distance between the Galactic center and the solar system, and r 
is the distance from the Galactic center. In studying the propagation of the cosmic-ray 
e^, we adopt the shape of the diffusion zone used in the so-called MED propagation 
model [36]; the diffusion zone is approximated by a cylinder with the half-height of 
L = 4 kpc and the radius of 7? = 20 kpc. Notice that the solution of the diffusion equation 
is insensitive to the diffusion coefficient Ke{E) for high energy [|3lll . as we will 
discuss in the following. In addition, the energy loss rate Z?ioss is given by the sum of 
the contributions from the synchrotron-radiation and the IC processes: b^^^^^^^\Ee,x) = 

^synch''^^(^e'^) + ^ic^'^''^'* ^) • For the Calculation of b'^^^^^\ we approximate that 
the strength of the magnetic flux density B is independent of position in our Galaxy; 
then we obtain 

b^^nr\Ee) = o^f,B\ (5) 

with 7e = Eel me and ax being the cross section of the Thomson scattering. In our 
numerical study, we use 5 = 3 /iG. Furthermore, b^^^^^"^ is given by 



^(Galaxy) (^^^^^^ ^ ^ J dEydEy^^ {Ey - Ey^^ ) ^^/tbg (^Tbg ' ^) ' (6) 



where c is the speed of light and the differential cross section for the IC process is 
expressed as [,37,1 



dOic 3 ax 



dEy AylE^ 



{TeqY{\-q) 

2(1 +r,^) 



(7) 



with Eg = AyeEy^^/nie, q = Ey/re{Ee —Ey), and fy^^ is the spectrum of the background 
radiation photon. Kinematically, 1 /Ay^ < ^ < 1 is allowed. The background photon in 
our Galaxy has three components: (i) star light concentrated in the Galaxy, (ii) star light 
re-scattered by dust, and (iii) the CMB radiation. The spectrum of the CMB radiation 
is isotropic and well known, while those of the first and second components depend on 
the position. We use the data of interstellar radiation field provided by the GALPROP 
collaboration [[bsI] . which is based on L39], to calculate /^^ in our Galaxy. 

The typical propagation length of electron per time scale of the energy loss is es- 
timated to be 0(0.1 kpc) for Eg ~ 100 GeV, and it becomes shorter as the energy in- 
creases. (Here, we have used the diffusion coefficient suggested in the MED propagation 
model, K{E) = 0.0112 kpc^/Myr x {Ee/l GeV)^-^^.) Then, the spectrum at the po- 
sition X is well approximated by 

%ss (^^'^) ^™mDM JEe dE'^ 

In our numerical analysis, we adopt the above approximated formula for the spectrum 
in our Galaxy. Then, y-ray flux from direction (b, I), where b and / are Galactic latitude 
and longitude, respectively, is obtained by line-of-sight (l.o.s) integral of y-ray energy 
density per unit time and unit energy as 

^f'''^''^\b,l) = ^[ dlLiciEyJ), (9) 



4k 



o.s 



where 



LiciEy.x) = c J dEedEy^J^fy^^{Ey^^,x)ff^''^''''^\Ee,x). (10) 



In studying the y-ray from cosmological distance, we need to understand the 
spectrum in the extra- Galactic region. In such a region, the spectrum is independent 
of the position, so the spectrum should obey 

djA^ = ^^^^^^^ + ^ [^ioss(^£.)/.(^£.)] + e(^£.), (11) 

where H is the expansion rate of the universe. Contrary to the case in our Galaxy, only 
the IC process with the CMB radiation contributes to the energy-loss process. Since the 
typical energy of the CMB radiation is so low that becomes non-relativistic in the 
center-of-mass energy of the IC process. Then, taking into account the red-shift of the 
CMB radiation, the energy loss rate is given by 

C^"^°^(^^.)=^^TT^pS(i+^r, (12) 

where P^mb^ — 0-26 eV/cm^ is the present energy density of the CMB, and z is the 
red- shift. 



The typical time scale of the energy loss due to the IC process is estimated to be 

^e/b[^°s"^°\ and is of the order of 10^"^ sec for Ee = 100 GeV (and becomes shorter as 
Eg increases). Because the energetic loses its energy via the IC process before the 
energy is red-shifted, we neglect the terms of 0(H) in Eq. (fTTI) and obtain 



1 



(now) 



.(Cosmo). „ 
^loss U.^e 



dE, 



,dNe 



dE 



(13) 



where P^m* — 1-2 x 10^^ GeV/cm^ is the present energy density of dark matter. Then, 
taking into account red-shift of scattered photon spectrum and the dilution due to the 
expansion of the universe, we obtain the y-ray from cosmological distance as 



where 



4> 



(Cosmo) 



C 

An 



dt 



1 



\Lic{t,Ey), 



(14) 



Licit, Ey) = {l+z)c J dEedEy^^ 



ACMB), „ w(Cosmo), p 



E'^{l+z)Ey 



(15) 



with fl^^^\t , Eyg^) being the spectrum of the CMB radiation at the time t. Notice that 
the astrophysical uncertainty is small in the extra- Galactic contribution, as is obvious 
from the above expression. 



NUMERICAL RESULTS 

Now, we are at the position to show our numerical results. First we will discussthe case 
where dark matter dominantly decays into /i^/i^ pair as an example. In such a case, 
energetic e"^ are produced by the decay of jj,^. Then, the positron fraction can be in good 
agreement with the PAMELA result while the total (e+ + e^) flux can be consistent with 
the Fermi data if moM ~ 0(1 TeV) and Tdm ~ 0(10^^ sec). In order to choose best-fit 
value with PAMELA data for Tdm in a given moM^ we calculate x^- In this anlysis, we 
simulate background flux of by the use of GALPROP code [3 SI . In the simulation. 



we adopt "conventional" model which is presented in [|4O0. 

In Fig. [B we show the simulated cosmic y-ray in the IC process. Here, we chose to 
take moM = 250 GeV, 500 GeV, 1 TeV, 2 TeV, and 4 TeV, and the lifetime is taken as 
Tdm = 1.0 X 10^7 sec, 6.2 x 10^^ sec, 3.4 x 10^^ sec, 1.7 x 10^^ sec, and 1.5 x 10^^ sec, 
respectively. For Galactic contribution, we plot the averaged flux mlO° <b < 20°. One 
can see Galactic contribution is much less than the observed IGL flux. On the other 
hand, the flux from extra- Galactic region makes a large contribution, especially when 
fnDM ^ 4 TeV, and it becomes comparable to the observed isotropic component of y- 
ray. Since the isotropic 7-ray flux observed by Fermi contains the extra- Galactic y-ray 
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FIGURE 1. IC 7-ray flux in the case where dark matter decays to /x+jU^ pair. The solid Une is 
the flux from the cosmological distance, which is isotropic, while the dashed ones are Galactic 
contributions. For each line, we take (moM, ^dm) = (250 GeV, 1 .0 x 10^^ sec), (500 GeV, 6.2 x 
10^6 sec), (1 TeV,3.4 x lO^^ sec), (2 TeV, 1.7 x 10^^ sec), (4 TeV,8.5 x 10^^ sec) from left to 
right. Here, we also plot observation data by EGRET ll4lll and Fermi ll33ll34ll . 
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FIGURE 2. IC 7-ray flux in the case where gravitino dark matter dominantly decays to 
second generation lepton. For each line, we take (moM, ^dm) = (250 GeV, 4.0 x 10^^ sec), 
(500 GeV,2.7 x 10^^ sec), (1 TeV,2.0 x 10^^ sec), (2 TeV, 1.8 x 10^^ sec), (4 TeV, 1.6 x 
10^^ sec) from left to right. 
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FIGURE 3. Primarily produced 7-ray flux in the case where gravitino dark matter dominantly 
decays to second generation lepton. The parameters are taken as the same as Fig. |2] 



and those from the other sources which are not identified, the extra-Galactic component 
should not exceed the observation; thus, such large dark matter mass may be constrained. 

For comparison, we also calculated IC-induced y-ray for scenario in which gravitino 
1//^ becomes dark matter. In the scenario, gravitino mainly decays to Wl, Zv, and 
hv. Then, is produced from charged lepton /, as well as weak- or higgs bosons; 
such can become a source of the cosmic ray [|42n to explain PAMELA anomaly 
[H, Si. Furthermore, total (e+ + e^) flux can be consistent with Fermi data. Even 
in this scenario, high energy y-ray is inevitably by e"^ in IC process. In Fig. [2l we 
show IC-induced y-ray in the scenario. Here, we especially consider the case where 
gravitino dominantly decays to second-generation lepton, and plot for moM = 250 GeV, 
500 GeV, 1 TeV, 2 TeV, and 4 TeV, by choosing best-fit lifetime with PAMELA as 
Tdm = 4.0 X 10^6 sec, 2.7 X 10^6 sec, 2.0 X 10^6 sec, 1.8 x lO^^ sec, and 1.6 x lO^^ sec, 
respectively. One can see that the Galactic and extra-Galactic fluxes from dark matter 
is much less than the present IGL- and isotropic-flux observations. This result is the 
consequence of softer spectrum from the decay of weak- or higgs bosons, compared 
to the previous /i^/i case. Besides, in i//^-DM scenario, the direct production of 
energetic photon by the dark matter decay is not negligible. In the scenario, 7-ray is 
produced directly from the decay of i//^, or following weak- and higgs bosons' decay; 
those photons also become the source of high energy cosmic 7-ray. In Fig.[3l we plot the 
primarily produced 7-ray. Here, we consider the same decay scenario and parameters 
as in Fig. [2l One can see that the Galactic and extra-Galactic fluxes from dark matter 
is much less than the present IGL- and isotropic-flux observations. Thus, with the 
observation data which is currently available, the primarily 7-ray is not constrained. 
However, if data for higher energy region is obtained, the scenario may be tested 
or constrained. For the other dark matter scenario, the spectrum of the cosmic 7-ray 



depends on the properties of the decaying dark matter. Detailed discussion for several 
specific dark matter models will be given elsewhere [44] . 

CONCLUSION AND DISCUSSION 

In this paper, we have studied the IC-induced y-ray, assuming that the PAMELA 
anomaly in the e+ fraction is due to the decay of dark matter. We have calculated the 
Galactic and extra- Galactic contributions separately, and shown that both of them are 
important. In our study, we have considered the case that dark matter dominantly de- 
cays into ji^ji^ pair, as well as the case of gravitino dark matter. However, as far as the 
PAMELA anomaly is explained in the decaying dark matter scenario, the IC-induced 7- 
ray flux is expected to be of the same order irrespective of dominant decay process. This 
is because, in order to explain the PAMELA anomaly, energetic should be produced 
by the decay, which induces the IC process. In our analysis, we found that large dark 
matter mass over nearly 4 TeV might be constrained by Fermi observation in the case 
the dark matter dominantly decays to /i^/i^ pair, while gravitino dark matter scenario 
is not constrained by the present observation. 

It should be also noticed that the production rate of energetic is extremely sup- 
pressed in the extra- Galactic region in another scenario of explaining the PAMELA 
anomaly, like the pulsar scenario. Furthermore, it is notable that the procedure to cal- 
culate this contribution has less uncertainty, and that the flux is isotropic. Therefore, if 
detailed data about the y-ray spectrum becomes available for various directions by the 
Fermi experiment, it will provide a significant information about the unstable dark mat- 
ter scenario. In particular, study of the y-ray flux from directions away from the Galactic 
center is important. For such directions, the 7-rays from the Galactic activities are sup- 
pressed. Thus, the background flux is expected to be comparable to or even smaller than 
the extra- Galactic contribution of the IC-induced 7-ray flux. 

At present, Fermi experiment provides the isotropic y-ray, by analyzing the flux data 
off Galactic center region [IssI]. According to their work, however, the isotropic compo- 
nent may contain contribution from unknown source in Galaxy, as well as extra- Galactic 
y-ray, and it is mentioned that more detailed analysis is needed in the future work. There- 
fore, if the y-ray spectrum for extra- Galactic component is precisely determined, we may 
see a signal of the IC-induced y-ray. Furthermore, it can be one of the possible tool for 
test of dark matter scenario. Also we note that the measurement of the flux from direc- 
tions close to the Galactic center, which is sensitive to the Galactic contribution, is also 
important. 
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